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MOLLIFICATION IN STRONGLY LIPSCHITZ DOMAINS WITH 
APPLICATION TO CONTINUOUS AND DISCRETE 
DE RHAM COMPLEXES^ 

ALEXANDRE ERN+ AND JEAN-LUC GUERMOND* * * § 

Abstract. We construct mollification operators in strongly Lipschitz domains that do not invoke 
non-trivial extensions, are L p stable for any real number pG [1, oo], and commute with the differential 
operators V, Vx, and V-. We also construct mollification operators satisfying boundary conditions 
and use them to characterize the kernel of traces related to the tangential and normal trace of vector 
fields. We use the mollification operators to build projection operators onto general H 1 -, H( curl)- 
and H (div)-conforming finite element spaces, with and without homogeneous boundary conditions. 
These operators commute with the differential operators V, Vx, and V-, are L p -stable, and have 
optimal approximation properties on smooth functions. 
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1. Introduction. Smoothing by mollification is an important tool for the anal¬ 
ysis and the approximation of partial differential equations. This tool has been in¬ 
troduced by Leray [H, p. 206], Soboleff 31, p. 487], and Friedrichs 18, p. 136-139]. 
The mollifying operation commutes with differential operators and converges opti¬ 
mally when the function to be smoothed is defined over the entire space R d . These 
properties may not be easy to achieve if the function in question is defined only in 
a non-smooth domain D , say only Lipschitz, since mollification by convolution as 
originally introduced in the above references requires to extend the function out¬ 
side D, which is a non-trivial task in general unless the boundary and the function 
are reasonably smooth. The boundary difficulty has been overcome by Blouza and 
Le Dret Q and Girault and Scott [lj| by redefining mollification using a convolution- 
translation technique so that mollification does not require information outside of D. 
The question is, however, more subtle when dealing with vectors fields where normal 
or tangent boundary conditions are involved, and, to our knowledge, has not yet been 
fully addressed in the literature. 

The first objective of this paper is to revisit the theory of mollification for scalar- 
and vector-valued fields in strongly Lipschitz domains with the following goals in 
mind: the mollification operators must be compatible with the De Rham complex 
(i.e., they must commute with the standard differential operators V, Vx, and V-), 
be L p -stable for any real number p £ [1, oo], and have strong convergence properties 
in the entire domain. Using a partition of unity subordinated to a covering of the 
boundary as done in Blouza and Le Dret @ and Girault and Scott [jj| is (seemingly) 
incompatible with the first constraint. The route that we propose instead is based on 
a shrinking technique of the domain using transversal vector fields. We also devise 
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a second sequence of mollification operators based on the extension by zero; the 
corresponding operators have the property that the mollified functions are in the 
kernel of the canonical trace operators in if 1 , iT(curl), and H( div) (i.e., zero trace, 
zero tangential trace, and zero normal trace, respectively). The main results of the 
first part of the paper are Theorems 13.31 and 14.41 and their respective corollaries. As 
an application of independent interest, we give in Theorem 14.71 a clear characterization 
of the kernel of the traces associated with the divergence and the curl operators. 

The second objective of this work is to use the mollification operators introduced 
in the first part of the paper to construct quasi-interpolation operators onto general 
finite element subspaces of H 1 , H( curl), and if (div), with and without homogeneous 
boundary conditions. We want these operators to satisfy three key properties: (i) 
Testability for any p £ [l,oo], (ii) commutation with differential operators, (iii) 
preservation of functions in the finite element spaces. Operators with such properties 
are important tools in finite element exterior calculus, see Arnold et al. a §5-4], 

§5.3], where they are termed bounded cochain projections. In particular, the 
above properties imply that the quasi-interpolation error is bounded by the best 
approximation error. 

The bases for constructing stable, commuting, and quasi-interpolation projec¬ 
tions have been laid out in Schoberl [28], [30| and Christiansen 0 , where stability 
and commutation are achieved by composing the canonical finite element interpola¬ 
tion operators with some mollification technique. Then, following Schoberl j2£j, the 
projection property over finite element spaces is obtained by composing these opera¬ 
tors with the inverse of their restriction to the said spaces. An important extension of 
this construction allowing the possibility of using shape-regular mesh sequences and 


boundary conditions has been achieved by Christiansen and Winther [12]. Further 


variants of this construction have recently been proposed. For instance in Christiansen 


11], the bounded cochain has the additional property of preserving polynomials lo¬ 


cally, up to a certain degree, and in Falk and Winther 0 it is defined locally. In the 
present work, we revisit the results of 0 by invoking our shrinking-based mollifica¬ 
tion operators which do not require extension outside the domain. We also present 
the results in the language of numerical analysis to make them accessible to a wide 
audience. The main result of this second part is Theorem 16.41 As an application, we 
give in Theorem 16.51 discrete Poincare inequalities for vector-valued functions. 

The paper is organized as follows. In ]J2] we introduce a shrinking technique of 
D that avoids difficulties with the boundary. We also introduce an expansion tech¬ 
nique. In ]J3] we use the shrinking technique to devise mollification operators that 
commute with differential operators, are stable in any L p space, and have approxi¬ 
mation properties on smooth functions. Using the expansion technique from ][2] we 
introduce in lj4] mollification operators that produce compactly supported functions 
and share the same properties as the shrinking-based operators. In we introduce 
the finite element setting that is necessary to construct canonical interpolation oper¬ 
ators on standard H 1 -, iT(curl)-, and iT(div)-conforming finite element spaces, with 
and without homogeneous boundary conditions. In [[5] we devise quasi-interpolant 
operators that are T p -stable, commute with differential operators, and preserve finite 
element spaces, with and without boundary conditions. 


2. Some geometry. In the entire paper, D is an open, bounded, strongly Lip- 
schitz, connected set in R d , int(D) denotes the interior of D , and D its closure. 
Points in R d and R d -valued functions and mappings are denoted using bold face; the 
Euclidean norm in M 9 , q > 1, is denoted |j.||^ 2 ( Rg ), or ||-||^2 when the context is unam- 
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biguous. We abuse the notation by using the same symbol for the ^-induced matrix 
norm and the norm of multilinear maps. 

2.1. Shrinking of strongly Lipschitz domains. Let D be a strongly Lipschitz 
domain, i.e., there are a > 0, f3 > 0, a finite number N of affine maps T n , n £ {1 :N}, 
and Lipschitz maps <h n : _B(0 R d-i, cc) — > R such that dD = Uti T n ({(x,z = 
& n (x)) | x £ H(0 R d-i,a)}), and for all n £ {l:iV}, 

T n ({(x,z) | x £ B( 0 R d-i,a), $„.(x) < 0 < $ n (x) +/3}) C int(D), 

T n ({(x, z) | x £ B(0 K d-i,a), 4> n (x) - /3 < z < $n(®)}) C M d \D, 

where S(0 R d-i,a) is the open ball of radius a in R d_1 centered at the origin. 

Since D is strongly Lipschitz and bounded, the combination of Theorem 2.7 and 
Lemma 2.2 from Hofmann et al. [22| implies that D has continuous globally transversal 
vector fields, i.e., there exist a vector field J £ C°(dD) and a real number 7 > 0 with 
the property that n(x)-J(x) > 7 at aspoint x on dD , where n is the unit normal vector 
pointing outward. Proposition 2.3 in [221 ] in turn implies the existence of a vector field 
j £ C°°(R <i ) whose restriction to dD is globally transversal and ||j( !E )ll <2 = 1 for all 
x £ dD. We then define the mapping: 

i/)j : l 1 * 3 3 : £ a: - Sj(x) £ R d . (2-1) 

Using Proposition 4.15 in [22j], together with the uniform cone property (see [2^, 
pp. 599-600]), we infer that there exists r > 0 such that 

ips{D) + B{0,6r)cD, VJ€[0,1]. (2.2) 


Lemma 2.1 (Properties of tps)- The following properties hold: 

(i) The map is of class C°° for all 6 £ [0,1]. 

(ii) For all i £ N, there is c such that ma|| D e (pg(x) — D l x ||^2 < c6, for all 
S £ [0,1], where D e denotes the Frechet derivative of order £. 

(iii) <ps(D) + B(0 , Sr) C D for all S £ (0,1]. 

(iv) The mapping x i-a x + t(ips(x) + ( Sry ) — a;) maps D into D for all t £ [0,1], all 
y £ B(0, 1) and all 5 £ [0,1]. 

Proof. The first two properties are consequences of j being of class C°° and D 
being bounded, while dm]) is just (12.21) . To prove (|rv|) . observe that t(<ps(x) — x) = 
<fts(x) — x. This implies that x + t(ips(x) + (Sry) — x) = ipts(x) + tSry £ <pts(D) + 
B(0,t5r) C D for all y £ B( 0,1), all t £ [0,1] and all S £ [0,1]. □ 

2.2. Expansion of strongly Lipschitz domains. Since D is bounded, there 
are xd £ R d and rn > 0 such that D C B(xd,td)- Let O = B(xD,rn)\D. The 
domain O is bounded, open, and strongly Lipschitz; hence, we can apply the above 
argument again, and deduce the existence of a vector field k £ C°°(R d ) that is 
globally transversal for 0, points outward D, and ||fc(a ;)|| / >2 = 1 for all x £ dO; note 
that dD C dO. We then define the mapping: 

'd 5 :R d 3xi—>x + 5k(x) £ R d . (2.3) 

As above, we infer that there exists £ > 0 such that 

i9 5 (0)+H(O,3<5C)c0, V<5 £ [0,1]. (2.4) 


Lemma 2.2 (Properties of ds). The following properties hold: 
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(i) The map fig is of class C°° for all S £ [0,1]. 

(ii) For all t £ N, there is c such that maxjg^ \\D tr dg(x) — D e x\\p < cS, for all 
S GjO,l]. 

(iii) 08 ( 0 ) + B(0,26() c O for all S £ (0,1]. 

Proof. The only novelty is the proof of (iii). Let x £ O, then there exists zg £ O 
such that (l+ 8 Lk)\\z 8 —x \\^2 < <5£, where L & denotes the Lipschitz constant of the field 
k in O. We observe that ‘dg{x) + B{ 0, 2 S() = 'dg(zg) + (# 5 (a;) — i?«(za)) + 13(0, 2<5£) C 
i9g(0) + B(0, 6 C) + B{ 0 , 25C) = 0s(O) + B( 0 , 3<SC) C O owing to □ 

3. Mollification without extension. We introduce in this section a mollifica¬ 
tion technique in strongly Lipschitz domains that does not require to invoke non-trivial 
extensions and that commutes with differential operators. The mapping ips : D —> D 
and r > 0 are defined in (ED and (12.21) . In what follows, Ja(a;) denotes the Jacobian 
matrix of ipg at x £ D. 

3.1. Mollification. Let us consider the following kernel 

p{y):=[ ^ if l|y||f2 < !’ (3.1) 

1 0, if H 2 /IU 2 > 1, 

where y is chosen so that f Rd p{y) dy = f B , Q ^ p(y) dy = 1. To be generic, we 
introduce q £ N, q > 1, and a smooth R 9X9 -valued field As : D —> R. 9X9 (related 
to the Jacobian JJa of ipg, see (13.51) below) such that for all l £ N, there is c such that 

sup ||.D ! (Ai(a;) — I )||£2 <cd, (3.2) 

x£D 

where I is the identity matrix in M qxq . Consider the following smoothing operators 
acting on f £ L 1 (U;R) and g = (gi,... ,g q ) T £ L 1 (D;R 9 ): 

(£f f)( x ) : = [ P{y)f{vs{x) + {8r)y)&y, (3.3a) 

Jb{ 0 , 1 ) 

(ICsg)(x) := Ag(x)(K.ggi(x),..., Kf s g q {x)) T . (3.3b) 

Property (lull) from Lemma 12.11 implies that 

ips{x) + ( 6r)y £ <ps(D) + (8r)B( 0, 1) C D, Mx £ D , My £ B( 0, 1). 

The means that the domains of /C® and ICg are indeed L 1 (U;R) and L X (D\\ R 9 ), i.e., 
there is no need to invoke extensions outside D. 

Lemma 3.1 (Smoothness). ICgg is in C°°(D;S. q ) for all g £ L 1 (D;'M. q ), and ICgg 
as well as all its derivatives admit a continuous extension to D. 

Proof. Owing to (13.21) and (I3.3bl) . and using the Leibniz product rule, it suffices 
to show that the statement holds for /C®. Let / £ L X {D). Let us prove first that /C®/ 
is continuous. Let x and 2 : be two points in D. Up to appropriate changes of variable 
we have 


ci/M - ACf/M = (dp J n (p(*^) - p) 1 ^)) /(y) dy, 


where we replaced <ps{x) + (8r)B( 0,1) and <fig(z) + (5r)B( 0,1) by D and used that 
p is zero outside the unit ball 13(0,1). The uniform Lipschitz continuity of p and tpg 








Mollification in Lipschitz domains 


5 


implies that there is c such that |p ^ y — p | < fp\\ x ~ z \\e 2 - As a 

result, we infer that 

\JC e s f(x) - ACf/(z)| < c{5r)- d - l \\f\\ L i (D )\\ x ~ z\\t 2 , 

which proves that AC®/ is Lipschitz continuous; hence AC®/ is uniformly continuous. 
This proves that AC®/ £ C°(D; K.) and AC®/ admits a continuous extension to D. Let 
us now evaluate the gradient of AC®/. Using the chain rule, we infer that 

V(A Cff)(x) = [ p{y)Sj{x)(Vf)(tp s (x) + (5r)y)dy 
Jb( 0 , 1 ) 

= jJ(*) / P(y)(yf)(ip s (x) + (6r)y)dy 
J B{ 0,1) 

= jJ(x)(5r) _1 f p(y)X7(f{cp s ( x ) + {6r)y))dy 
J B( 0,1) 

= f S/p{y)f(cp 5 (x) + (6r)y)dy. 

Jb( o,i) 

We can then conclude that V(AC®/) is Lipschitz continuous by using the same argu¬ 
ment as above and continue the argument by induction. 0 

3.2. Examples. Let / £ and g £ L 1 (Z?;M d ). Following (28| and [i(i 

§3], we define the following families of mollification operators: 


(£f/)(*) = 

/ p(y)f( t Ps{x) + ( Sr)y) dy, 

ffl(O.l) 

(3.4a) 

{lC 5 g){x) = 

/ p(y)^J( x )g{‘Ps{x) + (6r)y)dy, 

JB( 0,1) 

(3.4b) 

(*&)(») = 

[ p{y) det(I s (x))Sj 1 (x)g(ip 5 (x) + (Sr)y) dy, 

Jb( o,i) 

(3.4c) 

(*$/)(*) = 

[ p(y) det (S s ( x ))f(tp s (x) + ( Sr)y ) dy, 

ffl(O.l) 

(3.4d) 


for all x £ D. The superscripts in (13.41) refer to the fact that these operators are 
used to build projections onto finite element spaces that are conforming in the graph 
space of the gradient, curl, or divergence operator, or onto a broken finite element 
space (with no conformity requirement), see Theorem 16.41 below. The transforma¬ 
tions involving Jb are related to the classical Piola transformations. Furthermore, the 
functions lC c s g, IC d g, AC)?/ are of the form (I3.3bl) with A®(a;) = 1 and 

A^(a:) = lj (a;), A£(x) = det(JJi(a:))JI/ 1 (a:), A£(a;) = det(JJi(a:)) (3.5) 

Property (iii) from Lemma l2Tl implies that (13.21) holds true in the above three cases. 
Let p £ [1, oo]. Assuming d = 3, we define 


Z S ’ P (D) = {/ G L P (D) | 

V/ e L p (D)}, 

(3.6a) 

Z c - p (D) = {v£ L p (D) 

| Vxt; e L p (D)}, 

(3.6b) 

Z d ' p (D) = {v £ L P (D) 

| v-w £ L P (D)}. 

(3.6c) 


Lemma 3.2 (Commuting with differential operators). The following holds: 
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(i) V/Cf/ = ICgVf, for all f € Z**{D), 

(ii) Vx/Cfg = ACfVxg, for all g G Z C ’ P (D), 

(iii) V-lCjg = JCgV-g, for all g G Z d ’ p (D), 
i.e., the following diagram commutes: 


Z e ’ p (D) 

/Cf 

C°°{D) 


Z C,P (D) 

/q 

(D) 


Vx 


Z d ’ p (D) 


V- 


L p {D) 


Vx 


IC d 


C°°(D) 


V- 


)C h 5 


C°°(D) 


(3.7) 


Proof Upon setting T(x) = ipg(x) + 5ry, these identities are simple consequences 
of the chain rule: 


V(/oT)(*)=lJ(®)(V/)(T(*)), 

Vx(lJ (x)(goT))(x) = det(JJ^(a;))J^ 1 (a:)(V xg)(T(x)), 
Videt(I s (x))Sj\goT))(x) = det(J 4 (*))(V- ff )(T(*)). 

This completes the proof. □ 

3.3. Convergence. We now show that the smoothing operators defined above 
have interesting approximation properties. Owing to Lernma l2.il Jb and l^ 1 converge 
uniformly to the identity, and det(J,s) converges uniformly to 1. As a result, there 
is Jo £ (0)1] such that ||Ja — I||^a < |, ||JJ ^ 1 ||^2 < 2, and | det(JJq 1 )| < 2 d , for all 
J G [0, Jo] and all x G D. 

Theorem 3.3 (Convergence). The sequence (/C^rno^i is uniformly bounded in 
C{L P ]L P ) := £(L p (D-,M. q )-,L p (D;W 1 )) for all p G [0,oo]. Moreover, for p G [l,oo), 
\\l^sf — f\\LP(D;Ri) —■ > 0 as J —* 0 for all f G L P (D ; M 9 ). 

Proof. Owing to (13.21) and (I3.3bl) . it suffices to show that the statement holds for 
/Cf. (1) We show first that 1C f is uniformly bounded in £(L P ; L p ) by using the Riesz- 
Thorin interpolation theorem. The statement is evident for p = 00 with constant 
c = 1. Now consider / G L 1 (I?;R), then 


||/Cf/|| L i (D) < [ [ p(y)\f{ip s (x) + {6r)y)\dydx 
JD Jb( 0,1) 

< / p(y) / \f{<ps{x) + {8r)y)\dxdy 

JB{ 0,1) JD 

<[ p(y) f Ifi^Wdetilsiz))^ 1 dzdy < c||/|| L i (D) , 

JB( 0,1) J tp s {D)+(5r)y 

since J < Jo- The Riesz-Thorin interpolation theorem implies that \\lCgf\\LP(D) — 

C'||/||iP(D), so that \\m\\Lp(D) < ci||/||iP(D), with Cl = max(l,c). 

(2) Assume first that / is smooth over D , say uniformly Lipschitz with Lipschitz 
constant Lf, i.e., | f{x) — f(z)\ < Lf \\x — z\\ e 2 . Then, 


|£f/(*) - f( x )\ = 


>B( 0 , 1 ) 


P(y)(f(¥>s(x) + ( Sr)y) - f(x)) dy 


< / p{y)L f \\ips(x)-x + (Sr)y)\\ £ 2 dy <cL f S. 
Jb( 0 , 1 ) 
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In conclusion, there is Co = cmax(l, |Z?|) such that ||/Cf/ — J\\lp(d) < CqLjS. 

(3) We conclude by using a density argument and the triangle inequality. Let / £ 
L P {D). The space of uniformly Lipschitz functions being dense in L P {D ), there is a 
sequence of uniformly Lipschitz functions (/ n )ne n such that \\f n — /||lp(d) —> 0 as 
n —> oo. Then 

\\IClf - f\\ LP (D) < II £f(/ - f n )\\L*(D) + II ACf/„ - /„|| L p (13 ) + II fn - /||lp(D) 

< Cl II / - fn ||l p (L>) + CoLf^S + || fn - /||lp(D)- 

Let e > 0 and let n(e) be large enough so that ||/„( £ ) — /||lp(_d) < e. Setting So(e) = 
e/L/ n(e) , we have ||/Cf/ — /||lp(d) < e(ci +co + 1) for all 6 < 6o(e). In conclusion, for 
all e > 0, there is <5o(e) such that ||/Cf/ — f\\LP{D) < e(ci + Co + 1) for all S < <5 0 (e), 
which proves that ||/Cf/ — /||lp(d) —> 0 as S —> 0. □ 

We now can state a result that shows that the above smoothing technique is su¬ 
perior to mollification alone, i.e., contrary to the result originally stated by Friedrichs 
(see e.g., Thm 9.2]), strong convergence on the derivatives now occurs over the 
entire domain D. 

Corollary 3.4 (Convergence of derivatives). Let p £ [ l,oo). Then, 

lim || V(/Cf/ - /)\\ LP{D) =0, V/ £ Z**(D), (3.8a) 

O —^(J 

Hm||Vx(^- fl )|| i p (D) =0, \/geZ c ’ p (D), (3.8b) 

5—>0 

lim || V-(/Cfo - g)\\ L p(D) =0, V fl £ Z d ’ p (D). (3.8c) 

< 5—>0 

Proof. Using Lemma T3. 2 1 we infer that V/Cf/ = /CfV/, and Theorem 13.31 implies 
that /CrVf —> V f in L P (D ) as S —>■ 0, which proves (13. Sail . A similar argument holds 
for (liObl) and flUB . □ 

3.4. Convergence rate. We now establish convergence rates. 

Theorem 3.5 (Convergence rate). There is c such that 


\\fcsf ~ /||lp(D;R9) < Ci5 S |/|w'>.p(D;R9), (3-9) 

for all f £ W S,P (D; R 9 ), all S £ [0,<5o], and all s £ (0,1], p £ [l,oo) or s = 1, 
pe[i,oo]. 

Proof. Owing to (13.21) and (13.3bl) . it suffices to show that the statement holds for 
/Cf. Assume first that p < oo. 

(1) Let / £ W S,P (D) with s £ (0,1). We estimate /Cf / — / in L P (D) as follows: 


W^f-f\\ P L H D) = 


r D 


'£( 0 , 1 ) 


p(y)(f(<Ps(x) + ( Sr)y) - f(x)) d y 


dx 


< c 


J J lsp+d \\Vs(x) + (6r)y-x\\^ +d dxdy. 

Jb( 0 ,i)Jd \\ips(x) + {dr)y - x\\ e % 


Let us make the change of variables i?(0,1) 3 y z = tps(x) + ( Sr)y £ ips(D) + 
SrB(0 ,1) C D. Observe that the Jacobian of this transformation is bounded from 
above by Sr and 


II <Ps(x) + ( 5r)y - x\\p < \\<ps(x) - x \\ e 2 + Sr\\y \\ e 2 <cS. 
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Hence, 

II'CJ/-/e.(D, < c r p + d s- d J D JJj ( jz^£rr dxiz s 

(2) Let / € VL 1 ,P (-D). By proceeding as above we infer that 

]|£f/- /IIIp(d) ^ c / / !/(¥’«(*)+ (^)y) -/(*)| P dardy. 

JB(0,1) Jd 

Let us fix y G 5(0,1) and define the mapping xps '■ D 3 x i-a <pg(x) + (<5r)y G 
<ps(D) + 8rB(0, 1) C D. Observe that 

\\xfs(x)-x\\ P < \\v 5 {x)-x\\p+8r\\y\\p < cS , \\Dip s (x)-I\\ = ||.Dv? ( j(®)-I|| < cS, 

and x + t(if)s(x) — x) = x + t(tps(x) + dry — x) G D , i.e., 1/5 satisfies the assumptions 
of Lemma 13. 71 below. Hence, f D \ f(tp$(x) + (8r)y) — f(x) \ p dx < c8 p \\S7f\\ p LP ^ D y We 
conclude that ||/Cf/ - f\\ L p(D) < c6\\V f\\ LP{D )- 

(3) The case s = 1, p = 00 is treated similarly to (2). □ 

Corollary 3.6 (Convergence rate on derivatives). Let s G (0,1), p G [l,oo) or 
s = 1, p G [1, 00 ]. Then, there is c such that 

||V(/Cf/ - /)\\ LP{D) < cS s \\/f\ W s, P{D) , V/ G {x G L P (D) | Vo G W S ’ P (D)}, 

\\Vx(JC c s g - g)\\ LP{ D) <cS s \Vxg\ W s, P(D ), Vg G {v G L P (D) \ Vxo G W S ’ P (D)}, 

living-g)\\ L *(D) <cS s \X7-g\ W s, P(D) , Wg G {v G T p (5) | V-« G W S ' P (D)}. 


Proof. Let / G {0 G L P (D) \ Vo G VP S ’ P (-D)}, then 

l|V(/Cf/ - /)|| iP(c) - \\K%Vf - Vf\\ LP(D) since V/Cf = /QV 

< C(5 s ||V/||vi/s,p(d) owing to Theorem 13.51 

Proceed similarly for the two other estimates. □ 

Lemma 3.7 (Approximation). Let Ao > 0, and assume that : D — > D is a 
diffeomorphism of class C 1 such that \\ipx(x) — xjf^ < c' A and || D'if\{x) —1|| ^2 < A 
/or all x G D and all A G [0, Ao] - Assume also that the mapping g,x,t : x >->■ x + 
t(ip\(x) — x) maps D into D for all t G [0,1] and all A G [0, Ao]. Then, there is c 
such that the following holds: 


11 / f\\LP(D) < cA||V/||i,p(d), 

for all A G [ 0 , Ao], all f G W 1,P {D), and all p G [ 1 , 00]. 

Proof. (1) Assume first that / is smooth. Let x G D and v(t) := f(nx,t{x)) 
with t G [ 0 , 1]. The chain rule implies that v’(t) = Df(fj,x jt (x))(if}\(x) — x), thereby 
showing that 

f(ifx{x)) - f(x) = f v’(t)dt = f Df(nx,t(x))(il>x(x) - x)dt. 

Jo Jo 

Then, assuming that p < 00, we infer that 

II/O-0A - f\\ P LP{D) < ||^a(*) - x\\ p 2 ||V/(/x A)t (x ))|| p 2 dtdx 

<c'X p [ j \\Vf(g,\'t(x))\\ p 2 dtdx. 

Jo Jd 
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The assumptions on imply that the map is invertible and ||Z?/x A J ||^2 < 2, 
|det(D/x^()| < 2 d . As a result, 

II/o^a - f\\ P L P(D) - c ' XP f Q ^ ll V /(^)ll?2| de t(iA/x~^)| ci^dit, 

which finally implies that there is Co so that 11 / o xp x — f \ \ lp (D) < c o A 11 V/ 11 LP ^. The 
case p = oo is treated similarly. 

(2) If / is not smooth, we deduce from Corollary 13.41 that there exists a sequence of 
smooth functions converging to / in W 1,P (D), i.e., for all e > 0, there is a smooth 
function f t such that ||/ - / e || iP (D) < e and \\V f e \\ L p(D) < 2||V/|| iP(D) . Then 

11/ 01 h- f\\LP(D) < ll(/ — fe) O -lp\\\LP(D) + \\fe ° /e||z,P(£>) + II fe ~ /||lp(D) 

< ce + 2coA||V/|| lp(d) + e - 

The conclusion follows readily since e is arbitrary. □ 

4. Mollification with extension by zero. Note that the function ICsf defined 
in (|3.3bll does not satisfy any particular boundary condition. For instance, even if / 
is zero on dD , (ICsf)\go is not necessarily zero. Since preserving boundary conditions 
may be useful in some applications, we now construct a mollifier that has this property. 
Let C'g 0 (D;R <? ) denote the space of Revalued functions that are of class C°° and of 
compact support in D. Consider the mapping ds and the constant ( defined in (12.311 
and (|2.4D . Let IQ(a:) denote the Jacobian matrix of •ds at x £ D. 

4.1. Mollification. For any g £ L 1 (U;R 9 ), q £ N with q > 1, we denote by g 
the extension by zero of g over R d , i.e.,c(a:) = g(x) if x £ D and g(x) = 0 otherwise. 
Taking inspiration from Bonito et al. [7J, we introduce 

(IC e sof)( x ) : = [ P(y)f(&s(x) + (SQy)dy, (4.1) 

Jb( o,i) 

(K,s fi g)(x) := B 5 (®)(/Cf >0 p 1 (®), • ■ ■, K? so g q (x )) T , (4.2) 

for all x £ D, all / £ L 1 (L ) ; R), and all g = (g ±,..., g q ) T £ L X (D\ R 9 ), where B d - is a 
smooth R 9X9 -valued field (related to the Jacobian Kj of ds) such that for all l £ N, 
there is c such that 

sup ||l/(l> ( 5 (£c) — I )||^2 < cS. (4.3) 

x GD 


Lemma 4.1 (Smoothness and boundary condition). /Csp(g) is in Cq°(D-,W i ) for 
all g £ L 1 (D\ R 9 ) and all 5 £ (0,1]. 

Proof. The smoothness has already been proved in Lemma 13.11 Let k be the 
Lipschitz constant of the field k over D. Let es = (5£/(l + 5 k). Let a; £ I) be such 
that dist(a;, clD) < es ■ Then, there exists a point 2 £ dD such that dist(a:,;z) < es, 
i.e., 

ds(x) + B( 0, <5C) = d s (z) + B( 0, 5() + d s (x) - d s (z) 

= d s (z) + B( 0, <50 + x - z + S(k(x) - k(z)) 

C d s (z) + B( 0, 5() + B( 0, es + <5 ne s ) 

= ds(z) + B( 0, 5d) + (1 + <5k)£(s) = ds(z) + B( 0,2(50 
C i9s(0) + B(0, 2(50 C O, 
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owing to Lemma [2. 2 I' m). This implies that &s( x ) + (8Qy C O for all y £ 11(0,1), 
so that (/Cf o (f))(x) = 0 since + (SQy) = 0 for all y £ B{ 0 ,1). Hence, the 

support of /Cf 0 is compact in D. The same conclusion applies to K$, o- □ 

4.2. Examples. Let us proceed as in 33.21 Let Bf(®) = Kj(a:), Bf(ic) = 
det(Ka(a;))IK ( (" 1 (a:), and Bf (a:) = det(Ka(a:)). Lemma [2.21 )1) implies that (14.31) holds 
for these choices of Bj. Let / £ L^DjM) and g £ L 1 (D;R d ). We define the following 
families of mollification operators: 


0Cf,o/)(*) = 

[ P{y)f(&s{x) + (6()y)dy, 

Jb( o.i) 

(4.4a) 

(K^oflOO) = 

[ p{y)^s(x)g(x9s(x) + {SQy)dy, 

JB( 0,1) 

(4.4b) 

(^toflOO) = 

[ p(y ) det(K 5 (a:))Kj 1 (®)g(t? 5 (®) + (<5C)y) dp, 

ds(o.i) 

(4.4c) 

(^,o/)(*) = 

[ p{y) det{K s (x))f(ti s (x) + (<5C)y) dp, 

.1.5(0,1) 

(4.4d) 

x € D. Let p £ [1, oo]. If d = 3, we define 



Z«^(14) = {/ £ L p (H) | V/ £ X p (R d )}, 

(4.5a) 


Z C ’ P (D) = {v£ L P (D) | Vxn £ L p (R d )}, 

(4.5b) 


Z d ’ p (D) = {v£ L p (D) | V-v £ L p (R d )}. 

(4.5c) 


Lemma 4.2 (Commuting extension and derivatives). The following holds: 

(i) V/= V/^/or all f £ Z e,p (D), 

(ii) V xg =J7xg, for all g £ Z C ’ P (D), 

(iii) V -«7 = V g, for all g £ Z d ’ p (D). 

Proof Let / £ Z 6 ' P (D) and let xp £ C(j XD (K d ) be a (vector-valued) smooth function 
compactly supported in int(K d \H). Then, 


IR d 


xp-Vfdx = — / fV-xpdx = — / /V-0da: = O. 


ID 


Since xp is arbitrary, this proves that V/ is zero in R d \D. Now let xp £ Cff(D), then 


id 


[ xp-Vfdx = f /V-0 dx= [ , 
JR d Ifl 


— / xp-Vfdx = — / xp-Vfi lx = / fV-xpdx = / fV-xpdx = — / xp-Vfdx. 


id 


Since xp is arbitrary, this proves that (V/)m = V/. We have thus proved that 
V/ = V/. The argument for the other two equalities is identical. □ 

Lemma 4.3 (Commuting with differential operators). The following holds: 

(i) V/Cf >0 / = /q 0 V/, for all f £ Z**(D), 

(ii) Vx/q o g = /Cf 0 V xg, for all g £ Z^ P {D), 

(iii) V-/Cf i0 g = /Cf, 0 V-g, for all g £ Z d *(D), 
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i.e., the following diagram commutes: 


Z S ’ P (D) -——- Z C ’ P (D) -—-*- Z d ’ p (D ) -—-- L p {D) 

]jrc \r b 

^< 5,0 _ < 5,0 ^ < 5,0 ^< 5 , 

Cq°(D) -* CfT(D) — -» C£°(D) - CfT(D) 


(4.6) 


Proof. The proof is almost the same as that of Lemma T4.31 For instance, using 
the chain rule together with Lemma 14.21 we obtain 

V/Cf Q f(x) = [ p(y)K s (x) T Vf(#(x) + 6(y) d y 

JB( 0,1) 

= [ p(y)K s (x) T Vf(ti(x) + 6(y)dy = 1C c SQ X7f(x). 

Jb( o,i) 

Note here that it is critical that V/ = V/. The argument for the other two equalities 
is identical. □ 

4.3. Convergence. Similarly to H3.31 we can now state convergence results. 
Owing to Lemma [2.21 Ka and K^ 1 converge uniformly to the identity and det(Ka) 
converges uniformly to 1. As a result, there is <5o 6 (0,1] such that ||K^ — l\\p < i, 
||Kj 1 || <3 < 2, and |det(IK^ -1 )| < 2 d , for all S £ [0,<5o] and all x £ D. We combine the 
counterparts of Theorem 13.31 and Corollary 13.41 into one statement. 

Theorem 4.4 (Convergence). The sequence (/Ca,o), 5 g[o <5 0 ] uniformly bounded 
in C(L P ; L p ) for all p £ [1, oo]. Moreover, for all p £ [1, oo) 

Inn \\lC s ,of - /||lp( D; r«) =0, V/ G L P {D ; R«), (4.7) 

and 

lim ||V(/Cf 0 / - f)\\ LP(D) =0, V/ G Z^ p {D), (4.8a) 

o —>-U 

lim ||Vx(/C§ 0 g — g)\\LP(D) =0, Vg g Z c ’ p (D), (4.8b) 

< 5->0 

lim || V-(/Cf 0 g - g)\\ L , {D) =0, Vg € Z d ’ p (D). (4.8c) 


Proof. The proof of (14.71) is the same as that of Theorem 13.31 See the proof of 
Corollary [23] for the other three statements. □ 

Let s £ (0,1], p £ [l,oo) or s = 1, p £ [l,oo]. Let us denote by W s ’ p (D-,R q ), 
the space composed of the functions in W S ’ P (.D;R 9 ) whose extension by zero is in 
W s - p (R d ;R«). We set |/|^. lP(D . Rl) := \f\ W :*W)- 

Theorem 4.5 (Convergence rate). There is c such that 

W^Sflf - fWLP(D-Ri) < c $ S \f\w‘-r(D-,m)’ ( 4 - 9 ) 

for all f £ W s ’ p (D;R q ), all S £ [0,<5o], and all s £ (0,1], p £ ]l,oo) or s = 1, 
P € [l,oo]. 

Proof. The proof is identical to that of Theorem 13.51 □ 
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To state a convergence result using norms on D, we recall (see e.g., Grisvard [2(1 
Thm. 1.4.2.4, Cor. 1.4.4.5]) that 

W S ’ P {D ; IT) = Wq P (D; R 9 ) if sp ± 1, (4.10a) 

W s ' p (D;R q ) = W s ' p {D-,R q ) if sp € [0,1). (4.10b) 

(Recall also that the constants in the above norm equivalences depend on |sp — 1|.) 

COROLLARY 4.6 (Convergence rate on derivatives). Letp G [1, oo) and s G (0, |). 
Then, there is c (depending on \sp — 1\) such that 

l|V(/C® 0 / - f)\\ LP{D) < c6 s \\S7f\\ w *,p( D ), V/ G {r> G L P {D) \ Vv G W S ’ P {D)}, 

IIVx(/q. 0 fl - g)\\ L ,(D) < c6'\\Vxg\\ w ., r(D) , Vg€{ve L P {D) \ Vxv G W'*{D)}, 
l|V-(AC^ off - g)\\ LP (D) < cS a \\V-g\\ W :p(D), Vge{v€ L P {D) | V-« G W S,P {D)}. 


Proof. The proof relies on the commuting properties from Lemma 14.31 Theo¬ 
rem 031 and (14. lObl) . For instance, 

IIV(/C| 0 / - f)\\ LP{D) = || ICl 0 S/f - V /\\ LP(D) since V/C| 0 = /C^V 
<c<5 s |V/|^ a , V ^ IT) owing to Theorem 14.51 

< c' SiP 5 s \\Vf\\ W s, P{D) owing to (14. lObD . 

where c' s p depends on |sp — 1|. Proceed similarly for the two other estimates. 0 

Remark 4.1. The construction of 1C C S0 is similar in spirit to what has been pro¬ 
posed in Bonito et al. [?]. The curl estimates in Theorem 14.51 and Corollary 14.61 are 
identical to those in [3, Thm. 3.1]. 

Remark A.2. (sp > 1) Convergence rates on derivatives can also be derived for 
sp > 1, namely ||V(/Cf 0 / - f)\\ L p(D) < c5 s \S7 f\ W s, P ( D) for all / G L P {D) with 
V/ G W^ P {D), II Vx(/C§ 0 fif - g)\\ L p(D) < cS s \Vxg\ W s. P(D) for all g G L P {D) with 
Vxg G Wq’ p (D), and ||V-(/Cf 0 gr - g)\\ L p(D) < c5 s |V-g| ly »,P(D) for all g G L P {D ) 
with V-gG Wq’ p (D), where c depends on sp — 1|. Note that these estimates require 
boundary conditions on the derivatives. 

4.4. Traces of vector fields. In this section, we illustrate the use of the mollify¬ 
ing operator KLs.o- Let p G (1, oo). Recall the spaces Z C ’ P {D) and Z d ’ p {D) from (I3.6b[) - 
(I3.6cD . Since the trace operator 70 : W 1,p { D ) — )■ Wp’ p ( dD ) is surjective (see 
Brezis [s, p. 315], Grisvard 0, Thm. 1.5.1.2&1.5.1.6], McLean [ 25 ], Thm. 3.38] (for 
s G (4,|),p = 2)), letting {-,-)dD denote the duality pairing between W~p’ p {dD) 
and Wp’ p {dD), we define the bounded linear map 7 xn : Z C ’ P {D) —>■ W~p’ p {dD) by 

(7 xn{v),l)dD ■= / v-Vxw{l)dx - / w{l)-Vxvdx, (4-11) 

J D J D 

for all v G Z C ' P {D) and all l G Wp’ p {dD), where w{l) G W 1,p '{D) is such that 
7o(tc(i)) = l. Note that 7 X n( v ) = V |dDXn when v is smooth. The definition (14.111) 
is independent of the choice of w{l). Indeed, let W\,W 2 G W 1,p {D) be such that 
7o(tni) = 70(^2) = l , be., wi~w 2 G Wq’ p {D). Let (i/>„) nG n be a sequence in Cq°{D) 
converging to w\ — w 2 in Wq {D). Then, 0 = f D v-Vxi/>„ dx — f D xf n -Vxv dx, as 
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can be seen by replacing v by KT s v and passing to the limit 6 —> 0. Passing to 
the limit n —>- oo yields 0 = f D v-Vx{w 1 — w 2 ) dx — f D (w i — w 2 )-Vxvdx; hence, 
(7xn(v),jo(wi))dD = (7xnO)>7o02))aD, which establishes the claim. 

We also define 7 . n : Z d ' p (D ) —> W~p’ p {dD) by 

( 7 . n (v),l) dD := [ v-Vq{l)dx + [ q{l)V-vdx, (4.12) 

J D J D 

for all v £ Z d ’ p {D) and all l £ W~p' p ( dD ), where q{l) £ W 1,p (D) is such that 
70 (q(l)) = l, and {-,-}dD now denotes the duality pairing between W~p’ p {dD) and 
Wp ,p {dD). Reasoning as above, one can verify that this definition is independent of 
the choice of q(l). Note also that 7 - n {v) = V| QD' n when v is smooth. 

Let us now introduce 

Z c 0 p {D):=C™{D) ZC ’ P{D \ (4.13a) 

Z d ’ p {D) := Cg°{D) Z P(D) . (4.13b) 


Theorem 4.7 (Kernels of 7 xrl and 7 .„). Let p £ { l,oo). Then, 

Zq’ p {D) = ker( 7 XTl ), (4.14a) 

Z*’ P (D) = ker( 7 . n ). (4.14b) 

Proof. Let us do the proof for 7 Xn , the proof for 7 .^ is similar. 

(1) We first show that Z^ P {D) C ker( 7 xrl ), which is the easiest to establish. By 
definition there is a sequence of smooth functions (u n )raeN in C^°{D) converging to v 
in Z C ’ P {D). Let w be a function in C°°{D) (~l C°(D ), then 


0= / S7(wxv n )dx = / v n -S7xwdx — / ; 
J d Jd J D 


w-Vxv n dx. 


Both integrals on the right-hand side converge; hence, 


{lxn{v),7o(w))dD= / f-Vxiodx- / tu-Vxvdx = 0, 

Jd Jd 

for every function w in C°°{D ) n C°{D). This also implies that the equality holds for 
all w £ W 1,p (D), since C°°{D ) (7 C°{D) is dense in to £ W 1,p (D), see Lemma [■']. 1 1 
and Theorem 13.31 In conclusion, v £ ker( 7 Xn ) since 70 is surjective. 

(2) Let us now establish the converse, i.e., ker( 7 X7l ) C Z^ P {D). Let v £ ker( 7 XTl ). 
Since v £ Z c ’ p {D) C L 1 {D), v is differentiable in the distribution sense. Let t/> £ 
Cfi°(M d ), then 


(Vxr,^) = / D'Vxi/) dx = / v-'S/xif dx. 

J R d JD 

Using that v £ ker( 7 X7l ), the above equality implies that 

(Vx-u, t/>) = / r-Vxi/i dx = / xf-S/xv dx = / xf-S/xvdx. 


ID 
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This proves that Vxv = Vxv £ i 1 (R d ). Hence v £ Z C ' P {D). We can now apply 
(I4.8bl) from Theorem l4.4l since v e Z C ’ P (D ), i.e., the sequence (1C C S O i’) (5g [o j 0 ] converges 
to v in Z C ’ P (D). This proves that ker( 7 xrl ) C Z^’ P (D) since /C$ 0 v £ (see 

Lemma ED - □ 

Remark 4.3. (Zq = Z) The proof of Theorem 14.71 shows that Z^’ P (D) = Z C ’ P (D) 
and Z d ’ p (D) = Z d -P(£>); similarly, Zf p (D) = Z^ P (D). 

5. Finite element setting. We introduce in this section the finite element 
setting that we are going to use in the rest of the paper. We henceforth assume that 
D is a bounded polyhedron in R d . 

5.1. Meshes. Let ( Th)h>o be a shape-regular sequence of affine meshes. To 
avoid technical questions regarding hanging nodes, we also assume that the meshes 
cover D exactly and that they are matching, i.e., for all cells K,K’ £ Th such that 
K T K' and K fl K' T 0, the set I\ n K' is a common vertex, edge, or face of both K 
and K' (with obvious extensions in higher space dimensions). Given a mesh Th, the 
elements in K £ Th are closed sets in R d by convention. The following sets 

T k := {K’ £ % | K' n K ± 0}, (5.1a) 

Dk '■= int{a; £ D \ 3 K' £ Tk, x £ K'}, (5.1b) 

for all K £ Th, will be invoked in the following sections. The set Tk is the union 
of all the cells that touch I\, and Dk is the interior of the collection of the points 
composing the cells in Tk- 

We assume that there is a reference element K such that for any mesh Th and 
any cell K £ Th, there is a bijective affine mapping Tk ■ K —> K and an invertible 
matrix Ik € R dxd (not to be confused with JJ^) such that 

T K (x)-T K (y) =I K (x-y), Vx,y£K. (5.2) 

The shape-regularity assumption of the mesh sequence implies that there are uniform 
constants S, c such that 

|det(Jif)| = |A'||i?| -1 , Pa'II ti<<?h K , IUaII^ <c h h~ K \ (5.3) 

where Hk is the diameter of K. It can be shown that S' = — and S = —hp for 

Pk pk k 

meshes composed of simplices, where pk is the diameter of the largest ball that can 
be inscribed in K, is the diameter of K, and pg is the diameter of the largest ball 
that can be inscribed in K. 

5.2. Definition of S(x). In the arguments to follow, we are going to invoke 
smoothing operators like those defined in (J3l To avoid having to assume that the 
mesh sequence is quasi-uniform, we construct a meshsize function 1) £ C ILl (D: R) 
such that there are three uniform constants c, S , S' > 0 so that 

l|fl||wuoo(£, ;R ) < c, Sh K <h(x) <S'h K , Sx £ K, (5.4) 

for all K £ Th- The construction of this function is standard in the finite element 
literature. For instance, if the mesh is composed of simplices, consider the piecewise 
linear function whose value at any vertex of the mesh is the average of the mesh-sizes 
of the simplices sharing this vertex. 
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Following Christiansen and Whither [jjj], we introduce e £ (0,1) and define 

S(x) := ef)(a:), Vx £ D. (5.5) 

Then we can define ips and $5 like in (12.11) and (12.31) . and we can also define generic 
mollifying operators K.$ and lCs,o like in (I3.3bl) and (14.21) . Lemmas I2.1KJ2.2I hold for 
l £ {0,1} only, and the smoothness statement in Lemmas 13.1 FI4. 1 1 must be replaced 
by ICg(g) £ C' 1 (D;R 9 ) and K-s,o(g) £ Cq(.D;R 9 ) for all g £ L 1 (D;M. q ), respectively, 
since <5 is only Lipschitz. All the other statements in (J3]and fj4] remain unchanged. 

5.3. Reference and local finite elements. We are going to consider various 
approximation spaces based on the mesh sequence (7h)h>o and a fixed reference finite 
element (AT, P, S). We henceforth assume that P is composed of Revalued functions 
for some integer q > 1 and that P C W 1,00 (K\ R 9 ) (recall that P is a space of 
polynomial functions in general). The reference degrees of freedom and the associated 
reference shape functions are denoted {tri,..., <7„ 8h } and {$i,..., 0„ ah }, respectively. 
We denote N := {l:n s h} to alleviate the notation. We assume that the linear forms 
can be extended to £(P(AT);R), where V(K) is a Banach space such that 
V(K) C L l (K] R 9 ); see [THl, p. 39]. The interpolation operator Ig : V(K) —> P 
associated with the reference hnite element (AT, P, S) is defined by 

I^(v)(x) = ^2 &i{v)0i{x), Vx£K, Vv£V(K). (5.6) 

»e A/ - 


By construction, Xg. £ £(V(AT); P), and P is point-wise invariant by Ig. 

Let K be a cell in the mesh 7 h- We introduce a qxq invertible matrix Ak and 
define the mapping ipK £ £(L 1 (iF; R 9 ); L 1 (A'; R 9 )) by 


ip K {v) = A k (voT k ). 


ft can be shown (see 


15l . Prop. 1.61]) that upon setting 


Pk ■= {p = ^xiP) I £ e P}, 

'■= {&K,i}ieN S.t. <JK,i =&i ° 4>K, 


(5.7) 


(5.8a) 

(5.8b) 


the triple (AT, Pk,^k) is a finite element. Moreover, the interpolation operator 


l K (v)(x) = ^ u Kti (v)e K ,i(x), \/x £ AT, Vu £ V(K), (5.9) 

ieN 


where we have set Ok, i '■= i^ K l (^0- is such that Ik £ C(V(K)\Pk) and Pk is point- 
wise invariant by I K - Definition (15.8al) implies that Pk C W 1 ’ 00 (AT;R 9 ). More 
generally i/> K maps fW’ p (A:;R 9 ) to W l ’ p (K; R 9 ) for all l £ N, all p £ [1, oo] (with 
= 1, Vz > 0 if p = oo) and 


\' l / ;K \c(W , ’P(K;Ri);W l ’P(K;Ri)) — C ll-^fllt 2 l|Jic|l^ 2 |det(Jjf)| P, 
fe 1 ! C(W l -P(K-,Ri)-,W l -P(K-Ri)) — c ll^if lit 2 IU k 1 \\ 1 p |clet(JT jc)| p, 


for all K £ Th, (see e.g., 


13, Thm. 3.1.2] or 


15L Lemma 1.101]). 


(5.10a) 

(5.10b) 









16 


A. ERN, J.L. GUERMOND 


5.4. Structural assumptions. We henceforth assume that there is a uniform 
constant c such that 


II a a'|IH|Ar' 1 II£ 2 < c ||JJ*-||f 2 Pi < - 1 ||£2, (5-ii) 

so that, owing to (15.311 . 11A jc 11 £2 11A^- 1 11 ^2 is uniformly bounded with respect to K and 
h. We also assume that the degrees of freedom over K are either point values or 
integrals over edges, faces or K itself. This is formalized by assuming that 


|cj i (v)| < c 


ll«(2i)llf 2 (R9) 

IsFTI ll% 2 (R'0 ds, 


if point evaluation at 2j, 
otherwise, 


(5.12) 


where . is either an edge, a face, or K itself. All these mesh-related geometric 
entities are assumed to be closed sets. 

In the case of a point evaluation at 2^, we observe that since the cardinal number 
of £ is finite, there exists a distance Iq > 0 such that only one of the following 
situations occurs: (1) 2^ is a vertex of K; (2) 2^ is in the interior of an edge of K and 
is at least at distance Iq from any vertex; (3) 2^ is in the interior of a face of K and 
is at least at distance Iq from any edge; (4) 2 4 is in the interior of K and is at least 
at distance Iq from any face (with the obvious extension in higher space dimension). 

Let K £% and denote by {o,j}j £ M K the collection of points associated with the 
degrees of freedom in K defined by point evaluation. Note that there exists at £ K 
such that a,j = Tuidi) for all j £ A4 k- The shape-regularity of the mesh sequence 
implies that there is a constant c m i n (uniform with respect to j, K , and Th ) such 
that the open ball B(dj,c m i n hK ) has the following property: for all K' such that 
K' fl B(aj, c m \ n hK ) T 0 and every x £ K' fl B(dj, c m i n ha), the entire segment [x, dj] 
is in K'. An immediate consequence of this observation is that 


\\v(x)-v(dj )\\ e 2 < ||a:-aj||i 2 ||Vu|| L =o( J{ -/. R9 ), Vx £ K'f]B(dj, c m inh K ) ^ 0, (5.13) 

for all v £ Pk■ Note that this implies that B(dj,c ln i n hK) C Ta¬ 
in the rest of the paper, we define e max > 0 such that 


max max \\dj - {tpsfaAdj) + rS(dj)y)\\p < (5.14a) 

0 , 1 ) ^ 

uxeif (<?«(*)(*) + r5(x)B(0,1 )) C D K , (5.14b) 

for all K £ T h , all h > 0, and all functions S satisfying (15.51) for any e £ (0, e max ]. 

5.5. Finite element spaces. We introduce the broken finite element space 

P h (%) = {v h £ L^D-W) | ip K {v h \ K ) £ P , \/K £ %}, (5.15) 

where the statement ’•paiVh\ k) £ P is equivalent to Vh\K £ Pk- Notice also that 
P b (T h ) C W 1 '°°(7i;R ? ) := £ L°°(D\ R 9 ) | v\ K £ W 1 ’ 00 ^;®?), MK £ %} since 

Pk C W 1,00 (AT;R 9 ). We denote by : L P (D) -A P h (Th) the interpolation operator 
such that I b {v)\ K = Xk(v\k), f° r all K £ Th- 

We now introduce the notion of interfaces and jump across interfaces. We say 
that a subset F C D with a positive (d— l)-dimensional measure is an interface if 
there are distinct mesh cells Ki, K r £ Th such that F = dKi D dK r . We say that a 
subset F C D with positive (d— l)-dimensional measure is a boundary face if there 
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is a mesh cell AT G Th such that F = dl\ n dD. The unit normal vector rip on F is 
conventionally chosen to point from Ki to K r for an interface and to point outward 
for a boundary face. The interfaces are collected in the set T ^, the boundary faces are 
collected in the set T %, and we let Th = T^UT®. Let F G T° t be a mesh interface, and 
let Ki 7 K r be the two cells such that F = dKi n dK r ; the jump of v G W 1,1 (Th', R 9 ) 
across F is defined to be 

[v\p(x) = v\ Kl {x) - v\ Kr {x) aexGF. (5.16) 

Next we asume to have at hand a Banach space W C L 1 ( D ; R 9 ), with continuous 
embedding, where some notion of jump across interfaces makes sense. More precisely, 
we assume that there is a (bounded) linear trace operator jp ■ T / L 1 , 1 (iL; R 9 ) —> 
L 1 (9A";R t ), for some t > 1 and for all K € Th, and we define the notion of 7 -jump 
across interfaces as follows: 

Hf(®) = l Kl {v\ Kl )(x) - lK r {v\ Kr ){x) ae x G F. (5.17) 

We assume that |[u]J.(a:)| < |[u]i?(a;)|, for sea; € F, for all v G W l,1 (Th), so that 

v G R 9 ) =7 ([v]l = 0,VF6 T° h ). (5.18) 

We relate the notion of 7 -jump to the space W by assuming that 

vewnw 1 ’ 1 (%;&<!) =► ([«]?. = 0 , VF e T° h ), (5.19) 

and, conversely, that a function in W l '°°(Th', R ? ) with zero 7 -jumps across interfaces 
is in W. With this setting, we define 

P(Th) := P h (%) n W. (5.20) 

The above assumptions imply that 

P{T h ) = {v h € P h (%) | [v h ]l = 0, VF e T° h }. (5.21) 

Let F G T® be a boundary face and denote by Kp the unique cell such that 
F C dKp. We consider the global trace operator 7 : W 1 ’ 1 (U;M 9 ) —> L 1 (9U;R t ) 
such that 

i(v) ]f = 1Kf {v\k f ), VF&Tl (5.22) 

We assume that 7 can be extended to W into a bounded linear operator 7 : W —>• W d 
where W d is an appropriate Banach space, whose exact structure is not important 
for the time being. We define Wo = ker( 7 ), i.e., Wo = {v G W \ j(v) = 0}. Let us 
introduce P 0 {%) = P{%) n W 0 : 

Po(%) := {v h G P(Th) | l{v h ) = 0 }. (5.23) 

5.6. Examples. The present theory covers a large class of scalar- and vector¬ 
valued finite elements like Lagrange, Nedelec, and Raviart-Thomas finite elements. 
To remain general, we denote the three reference elements corresponding to the above 
three classes as follows: (AT, P g , S g ), (K 1 P C ,S C ) and (K, P d , E d ). The correspond¬ 
ing domains for the degrees of freedom are denoted V S (K), V C (K), V d (K). We think 
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of (P, P s , S s ) as a scalar-valued finite element (q = 1) and some of its degrees of free¬ 
dom require point evaluation, for instance ( K , P g , E s ) could be a Lagrange element. 
We assume that the finite element (P , P c , S c ) is vector-valued (q = d ) and some of its 
degrees of freedom require to evaluate integrals over edges. Typically, (P, P C ,P C ) is 
a Nedelec-type or edge element. Likewise, the finite element (P, P d , £ d ) is assumed 
to be vector-valued (q = d) and some of its degrees of freedom are assumed to require 
evaluation of integrals over faces. Typically, ( K,P d ,E d ) is a Raviart-Thomas-type 
element. The arguments developed herein do not require to know the exact structure 
of the above elements. 

The above assumptions imply that it is admissible to choose V s (P) = W S ’ P (K ) 
with s > |, V C (P) = W S ' P (K) with s > ^-, and V d (K) = W S ’ P (K) with s > | 
(recall that denoting by M a smooth manifold of dimension d! in P, the restriction 
operator to M is continuous from W S,P (K) to L P (M) provided s > ^y-)- Note that 
it is also legitimate to choose 

V s {K) = W d ’ 1 (K), V d (K) = W 11 {K), V c (K) = W d ~ 1 ’ 1 (K), 

since W (Ll (K) C C°(K)), functions in LF 1,1 (iL) have a trace in L l (dK ), and functions 
in W“ _1,1 (P) have integrable traces on the one-dimensional edges of P. 

Let ipx, xj) c K , xpft be the linear maps introduced in (15.71) for each of the reference 
finite elements defined above. In practice if )is the pullback by T K , and and ipp 


are the contravariant and covariant Piola transformations, respectively, i.e., 

Aff = 1, #)=»» T k , (5.24a) 

A c x = JT j ^-(«) = jT (t , D Ta . ); (5.24 b) 

A^- = det(lif) J^ 1 , = det(JTjf) o T K ). (5.24c) 

Note that c = 1 in (|5.11l) for the above examples. We consider the following y-traces: 

1k( v \k)(x) := V\ K (x), Vx G P, (5.25a) 

7 k( v \k){x) := v lK (x)xn F , \/x G P, (5.25b) 

1 k( v \k)(x) := K (x)-n F , Vx G P, (5.25c) 

and the following finite element spaces: 


P g (T h ) := {v h G L\D) | <P s K (v hlK ) G P g , VP G T h , K] g = 0, VP G P° h }, (5.26a) 
P c (T h ) := {v h G L\D ) | *) GP C , VP G %, = 0 , VP G 7£},(5-26b) 

P A (T h ) := {v h G L\D ) | eP d , VP G T ft , K]^ = 0, VP G J£}, (5.26c) 

where we simplified the notation by using {vhjp instead of etc. Note the 

conformity properties P s (T h ) C Z^ P (D), P c (T h ) C Z C ’ P (D ), and P d (7;) C Z d ’ p (D). 
Likewise, observing that Z g ’ p (D) := {v G Z S ’ P (D) \ j e (v) = 0} etc., we define 

P§(%)~ PZ(T h )nZl’ p (D), 

P§(%):= P c (%)nZ dp (D), 

P d (T h ) := P d (T h )nZ d ’ p (D). 


(5.27a) 

(5.27b) 

(5.27c) 
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We finally denote by X®, X d , X % and X® 0 , Z d 0 , X£ 0 the canonical interpolation 
operators associated with the finite element spaces P g (Th ), P c (Th), P A (Th), and 
Poijh), P§(%), Note that X® : W S ’ P {D) —► Ps{T h ) C L P (D) is stable 

provided s > ( j, Zf t : W S ’ P (D) —> P c (T h ) C L P (D) is stable provided s > 

(= | for d = 3), and Z £ : W S,P (D) — > P d {Th) C L P {D) is stable provided s > 
d ~ < 'p~ 1 ' > = p- We finally assume that the polynomial degrees in each of these spaces 
are compatible so that the following commuting properties hold with s > j } : 


V Vx , V- 

V s (D) - - V C {D) - ► V d {D) - * L p (D) 

XI Z c h Z d Z\ 

P g (Th) --—► P c (T h ) — P d {T h ) - P b (T h ) 

where 

V g {D) = {/ G W S ’ P (D) | V/ € W s ~p’ p (D)}, 

V C (D) ={j£ W s ~p’ p {D) | V xg G W s -p' p (D)}, 
V d {D) = {g€ W s ~p’ p (D) | S7 g G W s ~p’ p (D)}, 


(5.28) 


(5.29a) 

(5.29b) 

(5.29c) 


and Z b is an interpolation operator only involving integrals over mesh cells. Likewise, 
upon introducing 


V 0 S (D) = {/ G V s (D) | f ldD = 0}, (5.30a) 

V 0 C (D) = {g G V C (D) | gxn ]dD = 0}, (5.30b) 

V d (D) ={jG V d (D) | g-n ldD = 0}, (5.30c) 

we assume that the following diagram commutes: 

V 0 S (D) - V 0 C (D ) - V d (D) - L p (D) 

x m 

Pq(Th) - PS 

6. Stable, commuting, quasi-interpolation projection. We introduce in 
this section a family of finite-element-based quasi-interpolation operators (with and 
without boundary conditions) that are X p -stable, commute with the standard differ¬ 
ential operators V, Vx, and V-, and preserve the above finite element spaces. 

6.1. The operator ZhlCg. Owing to the properties of the smoothing operators 
established above, it makes sense to consider the discrete functions X®/C®/ G P s (Th), 
W e PS(T h ), Z^f G P h (T h ), zycig G P c (%), Xfaiqg G PS(%), Z d lC d 5 g G 
P d (Th) and Z d 0 K, d g G P d (Th) for any integrablc scalar-valued function / and any 
integrable vector-valued function g. We now establish some stability properties of the 
restrictions of the operators Xj)/Cf, 2f n /Cf, X!)/Cj(, X?/C?, Xf n /C?, X^/Cf, and X? n /rf to 
the discrete spaces P*{%), P g (%), P h {T h ), P c (T h ), P§(T h ), P d (%), and P$(T h ). 

To avoid repeating proofs seven times, we denote by Zh one of the interpolation 
operators introduced above and K,$ the corresponding smoothing operator; likewise, 
the range of Zh is denoted P{Th)- We assume that P(Th) is composed of Revalued 
fields. 


-'hO 


{%) 


Vx 


T d 

x h0 


T b 
x h 


PSKT h ) 


P h (Th) 


(5.31) 
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Remark 6.1. (Boundary conditions) Note that we do not invoke /Cf 0 , /C]] 0 , and 
ICg 0 in the above construction. The theory to be exposed in the next section holds 
by using K, f, 1C C S , and 1C) in all the cases, whether homogeneous boundary conditions 
are enforced or not in the discrete spaces. 

6.2. Testability of the operator IhICs. We start with a key result in the 
spirit of [HI, Lem. 4.2], see also [2fil . Lem. 6]. This result is crucial to devise a quasi¬ 
interpolation operator that preserves the finite element space P(Th)- 

Lemma 6.1 (Discrete T p -approximation). There is c sta b > 0, uniform with re¬ 
spect to the mesh sequence, such that || f h -Th>Csfh\\LP(D-Ri) < c sta .b^\\fh\\LP(D-Ri) for 
all e £ (0, e max ], all f h £ P(%) and all p £ [1, oo]. 

Proof. (1) Let f h £ P(%) and let us set e h := f h - l h ICsfh and e \=J h - /C s f h \ 
note that = fZ^e. Let K be a cell in Th , then using that 9k, i '■= we have 

pK,i\\LP(K-Ri) < det(l/f)p ||A^ 1 ||^2||0 i || iP( ^ ;]R9) for all i £ {l:n sh }, and we infer that 

\\ e h\\LP(K-,Ri) = W^hO\\ lp(K-, Ri) < ^ \ a K,i (e) | \\9k,i || LP(K-Ri) 

iGAf 

< det(JJ if )p||A ^ 1 ||^2 ^2 l <J A,i(e)|, (6.1) 

ieM 


The rest of the proof consists of estimating crx,i(e). 

(2) Let us assume first that the degree of freedom gk,i is a value at a point 
a,j := Tft-(ai) in K. Then using the assumption (15.121) and the definition (15.71) of 
i( k, we infer that |crR-,i(e)| < c||Ax||^||e(Oj)||^2. By proceeding as in the proof of 
Theorem 13.51 (step (2)), we obtain 

eK) = h(aj) - ICsfhiaj) = / p(y)(f h (a j ) - f h (<Rs( aj ){aj) + S(aj)ry)) d y. 

JB( 0 , 1 ) 

Owing to (15.131) and (I5.14al) (recall that e < e max ), we have 

ll e ( a l)llt 2 < c max || f h {aj) - f h {y>8( a . j ){ a j) + K a ])ry)\\p 
yeB(0,l) 

< c'd(a J ) max || Vfh\\L°°(K'-,m) < c"eh K max || Vfh\\L°°(K'-,m)- 
K'e/K K-'^Tk 

Finally using a local inverse inequality, which is legitimate since the mesh sequence 
is shape-regular, we infer that |crK-,i(e)| < ce||A^c||^ 2 \\fh\\L°°(D K ;R q )- Note that the 
purpose of the above argument is to account for the fact that fh is (a priori) only 
piecewise Lipschitz (i.e., can be discontinuous across interfaces) but fh is necessarily 
continuous at a,j. 

(3) If the degree of freedom UK,i is an integral over an edge, face or over K, we 
use (15.121) , i.e., Icr^^(e)| < c||A^|| f2 |^-j J Sk _ |je|| f 2 ds. We define Ts K<i = {K ' e 
Tk | Sk,i C K’} and we introduce S° Ki = {x £ Sx,i \ <Ps(x){x) + 5(x)rB(0,l) C 
7s k ,J and S^ ti = S K ,i\S° Ki . Then using (I5.14b|) and setting ip 5 (x,y) = ip s ^(x) + 
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5(x)ry, we have 

[ \\e\\pds< [ [ p(y)\\f h {ip s (x,y)) - f h (x)\\p dyds 

JS° K i •' S K, i •'#(0,1) 

/ eB( oi, \\fh{^s{x,y)) - fh{x)\\p dsdy 

z ' / y 6 B( 0 ,l) 

,i K , ^1 ~s k I ^8( x ’y) eK ' 

< c\SK,i\eh,K E n v ^iii “(_fsr';R9) < cis^ikll/hiu ~(D K ;R1), 

K 'eTs K ,i 

where we used the shape-regularity of the mesh sequence (i.e., hx' < cHk) and an 
inverse inequality. Note again that the above construction is meant to account for the 
fact that fh is (a priori) only piecewise Lipschitz. Moreover, if x £ Sj^ i , then there 
is y £ B(0 ,1) such that z := ip$( x )(x) + 5(x)ry is not in Ts Kii ; then mesh-regularity 
implies that c dist(x, dSie,i) < \\z — x^\jp and that \\z — x\\p < cS(x) < c'eliK- 
Combining these bounds, we obtain that \Sf c t \ < cehx\dSK,i\ < c'e\SK,i\ (with the 
convention that the O-dimensional measure of a point is 1). As a result, we infer that 

[ \\ e \\e 2 ds < f {\\fh\\t 2 + \\lCsfh\\e. 2 )ds 

J s e K ,i Js a Kti 

< c \\fh\\L°°(D K -,R*)\SK,i\ < c ' e l‘S'if,i|||// l ||L<»(D K ;R9)- 

Combining the above two estimates yields \crK,i(e.)\ < ce\\AK\\£ 2 \\fh\\L°°(D K ;S.i)- 

(4) We have established that | (JK,i{e)\ < ce|| KkWp\\ fh\\L°°(D K -Ri) for all possible 
degrees of freedom. Using the fact that || \\p || A^- 1 \\p is uniformly bounded together 

with an inverse inequality from L°°{DK' 1 ^- q ) to L p {Dk\ R 9 ), we deduce that 

\\fh - lhfcsfh\\LP(K-,m.Q) = l|e/i|Up(if ; R9) 

< cedet^i^)? HA^ 1 1|^21|Ax || e 2 \\fh\\L°°(D K -,Ri) 

< ce \\fh\\LP(D K -,RQ)- 

We infer the desired result by summing over K G Th and by invoking the shape- 
regularity of the mesh sequence. □ 

The above lemma implies that ||(I — IhK-5)\p(T h ) Wc(lp-lp) < c s t a be for all e £ 
(0, e max ]. From now on we choose e once and for all by setting e = e m i n with e m i n := 
min(e max , (2c s t a b) -1 )- Lemma f6.11 then implies that 

\\(I - l h ICs)\p^r h )\\c(LP-,LP) < 2' (6-2) 

This proves that IhK-s\p(j- h ) is invertible for this particular choice of e. Let Jh : 
P{Th) —> P{Th ) be the inverse of l h K. s |p(r h ), i-e., 

J h lhlCs\p(Th) = Phfc-6\P(T h )Jh = I- (6.3) 



Note that the definition of Jh implies that \\Jh\\c(LP-,LP) < 2. 

Lemma 6.2 (A p -stability). Let e = e m i n . There is c(e m i n ), uniform with respect 
to h, such that the following estimate holds: \\lThfcs\\c(LP;LP) < c(e m i n ). 
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Proof. Let / € L p (D;R q ) and assume p < oo. Then 


WWsfW 


V 

LP(D;E9) 



PK,i{K-5 f)0K,i{x) 


P 

P 


dx 


< c 


\ a K,z{1<:sf)\ P \\dK,i{x)\\ P e2 dx. 

K£T h Jk ie At 


Using (15. 1211 . we infer that 

IIWIII P ( W <c £ ^||Ak||^||/C,/||^( W ||A^||^|^|, 

KeT h iGAT 

since \<T K ,i(ICsf)\ < c\\A K \\p\\JC s f\\ L <*,( K . Rq) and \\O k A x )\\l^(k^i) < c||A^ 1 || £ 2 . We 
conclude by invoking Lemma 16.31 below. The argument for p = oo is similar. □ 

Lemma 6.3 (Local inverse inequality). Let e = e m i n . There is a uniform constant 
c > 0 such that 


ll^5/IU°°(Jf;R«) < Ce min |AT| p ||/||lp(Dk;R9)) 


(6.4) 


for all K £ %, all h > 0, and all f £ L P (D ; 1 

Proof. Let x £ K. Since the function p is bounded, we infer that 


\\ICsf(x )\\ e 2 <c f || f(<Ps( x )(x) + S(x)ry )\\ e 2 d y. 
Jb( o,i) 

The condition (I5.14bl) implies that 


\\IC S f(x)y<c\\S-TL^D K ) / ll/(*)ll^ dz < ce^ n hK a \D K \ L -p\\f\\ LnDKm , 

J Dk 

and we conclude using the shape-regularity of the mesh sequence. 0 
6.3. Main result. We now define the following operator 

Jh = JhlhlCs, (6.5) 


and we state the main result of this section. 

Theorem 6.4 (Properties of Jh). The following properties hold: 

(i) P(Th) is point-wise invariant under Jh; 

(ii) There is c, uniform with respect to h, such that \\Jh\\c(LP;LP) < c and 

11/ - Jhf\\LP(D;Ri) < C inf ||/ - /h ||lp(_D ; R9) , 
fh&P{Th) 


for all f £ L p (D;R q ); 

(iii) Jh commutes with the standard differential operators, i.e., the following diagrams 
are commutative: 


Z S,P (D) --—► Z cp (D) — -► Z d ’ p (D) -—-- L p {D ) 


7 s 

Jh, 

J c h 

& 

L, v 

L, v X 

v- 


V Vx , V- , 

P s (T h ) -- P c (%) -- P d (T h ) -- P h {T h ) 


( 6 . 6 ) 
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Z$’ P (D) 


Zq P {D) 


Vx 


Zo’ P (D) 


V- 


ji 

yjy 


hO 


Jl 


hO 


P$(T h ) 


P 0 c (Th ) 


Vx 


7 d 


PoiTh) 


L P {D) 


jh 


P h (T h ) 


(6.7) 


Proof. The first property is a consequence of Jh\P{T h ) = ^ since Jh°Jh = 
Jh\p(T h ) 0 Jh = Jh- The second property is proved by observing that the L p -operator- 
norm of Jh is bounded by 2 and that of Xh^s is also uniformly bounded, as established 
in Lemma since e is now a fixed real number. Moreover, using that Jhfh = fh 
for all fh G P(Th), we have 

11/- J4/||lp(D;Ri) = inf \\f ~ fh ~ Jh{f ~ fh)hp(D-m) 

fh&P(Th) 

< . inf ,(f + \\Jh\\c(LP-LP))\\f - fh ||l p (D;R 9 ) < C inf 11/ - fh\\LP(D;Ri), 
jh^r’ylh) jht^Uh) 

which establishes ©■ Let us now prove (|m)l . We are just going to show that the 
leftmost square commutes in the top diagram; the proof for the other squares is 
identical, and whether boundary conditions are imposed or not is immaterial in the 
argument. Let us first show that J^S/cfh = S7(J^(fh) for all (j>h G P s (7/). We observe 
that 

= V(2*/Cf |JJg(7i) Jjft/0 = V(Jf/Cf = 2£V(/Cf = l c h IC c 5 V(J^h), 

where we have used that I = J® (see (E3D) , then Ph^- 6 s\pi(r h )Jh ~ Ph^s^h 

(the range of J® is in P g (7/t)), followed by VZj) = I£V (see diagram (15.281) ') and 
V/C| = /QV (see diagram (13.71) ). Since V(J®</>/j) G P c (7/) (see diagram (15.281) 1. the 
above argument together with (16.31) proves that 

= (W |pc(Tfi) )v(j®^) = 

In conclusion, JfS/f’h = Now we finish the proof by using an arbitrary 

function cf> G V s (P). We have 

= j^v(/c^) = 4V(j®/cf/») = v(j®j®/c|</>). 

The last equality results from the fact that J^S/(j>h = V(J|</^) for all (fh G P®(7/), 
as established above. This proves that = VJffcj). □ 

Remark 6.2. (Approximation) Theorem 16. dl iil shows that the quasi-interpolation 
error is bounded by the best approximation error. Estimates of best approxima¬ 
tion errors in fractional-order Sobolev spaces have been obtained recently in Ern 
and Guermond [16] for general finite element spaces. As an illustration, consider a 
P c (7h )-based finite element approximation of a field A G Z C,P (D) (typically, with 
p = 2). Suppose that the natural stability norm for this problem is that of iL(curl) 
and that the finite element solution Ah G P c (Th) satisfies the a priori error estimate 
Hi*(curl) < cM ahGP c {Th) H-4. —a /j .|| Pr(curl) . Then, taking a h = JfA and using 
the commuting property leads to the bound 

11-4 — AhHj^curi) < c(|| A — l 7/A|| i 2( £3 ) + || Vx-4. — J' d (Vx J 4)|| £/ 2( £ ,)). 

Assume that A, VxA G H r (D) for some real number rG (0,H1] where k is the 
degree of the finite elements composing P c (7/). Then, using Theorem l6.4l ii') together 
with [H, Cor.5.4] leads to ||A - A h \\ jf( curl ) < ch r (\A\ H r( D) + IVxAl^^)). Note 
that no lower bound on r is assumed a priori. 
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6.4. Discrete Poincare inequalities. We illustrate the usefulness of the oper¬ 
ators constructed above by proving discrete Poincare inequalities for H (curl)-elements 
in dimension d = 2,3; we expose the material for d = 3. Assume also that D is par¬ 
titioned into M connected, strongly Lipschitz subdomains D i, • • - , Dm- We consider 
two piecewise-smooth second-order tensor fields c and qj, i.e., we assume that these 
fields are in 


W 1 ’ 0O (uf =1 A):={''eL CO (0) I V(i/, D J S [L°°(A)] d , * = (6.8) 

where L °°(E) := L°°(E- : M dxd ). We additionally assume that e and p are symmetric 
and the smallest eigenvalue of each of these two tensors is bounded away from zero 
from below uniformly over D. Consider the following Maxwell eigenvalue problems: 
Find E and 0 ^ to £ R. such that 

Vx(qj ~ 1 VxE)=uE, V-(c£J) = 0, Exn\ gD =0, (6.9) 

Find B and R such that 

Vx(qj -1 Vx.B) = ujB, V-(cB) = 0, (cB)-n\ gD = 0. (6.10) 

Upon setting H xrl := {z £ Ff(curl) | V-(£z) = 0, zxn^ gD = 0}, H n := {z £ 
H( curl) | V-(cz) = 0, {<s:z)-n\ gD = 0}, the L 2 -theory of the well-posedness of this 
problem is based on the following embedding inequality: There are c > 0 and s > 0 
(both depending on D and <e) such that 


||e||H.(r») < c||Vxe||i 2 ( D) , Ve £ H xn , (6.11a) 

l|b||jj»(D) < c ||Vx 6 || i 2 ( £3 ), Vb £ H. n , (6.11b) 


provided dD is connected and D is simply connected, respectively. The above in- 
equalities, proved in Bonito et al. [g], generalize classical inequalities established by 
Costabel 141 and Birman and Solomyak Q assuming that the tensor c is smooth over 
the entire domain. 

Let us consider the finite element approximation of the above eigenvalue problem 
using the setting described in the previous sections. The approximation theory of this 
problem is non-trivial, especially when using finite elements that do not fit the De 
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Rham Diagram. We refer to the book of Monk 
for an overview on the topic. 

Let P c (Th), Po(Th) be defined as above. A key step for approximating (16.91) 
or (16.101) consists of establishing the following discrete Poincare inequalities: There is 
c > 0 , uniform with respect to h, such that 


Il e /t||i, 2 (£)) < c ||Vxe/ 1 ,||x, 2 ( J c)), Ve/i £ Hh, xn , (6.12a) 

||b/i||.L 2 (D) < c ||Vxb h \\l 2 (d), V6/j £ Hh t - n , (6.12b) 

where H hiXn := {v h £ P§(Th) | f D (ev h )-S7q h dx = 0, Mq h £ P§(Th)} and H h . n := 
{vh £ P c (Th) I J D {ev h )-Vq h dx = 0, Mq h £ P g (%)}■ There are many ways of proving 
(16.12al) - (|6.12bl) when £ is smooth, since in this case it can be proved that the Sobolev 
index s in (16.11a!) is larger than tj . The first route described in [ 2 ]], §4.2] consists of 
invoking subtle regularity estimates from Amrouche et al. [2, Lemma 4.7]. The second 
one, which avoids invoking regularity estimates, is based on the so-called discrete 
compactness argument of Kikuchi [23j and further developed by Monk and Demkowicz 
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27] and Caorsi et al. [9|. The proof is not constructive and is based on an argument 


by contradiction. 

We now show that using the approximation operators described in the previous 
sections gives a direct answer to the above question without requiring any particular 
condition on the Sobolev index s in (16.11a!) : see also Arnold et al. [$, Thm 3.6]. 

Theorem 6.5 (Discrete Poincare). Assume that dD is connected (resp., D is 
simply connected). Then there is a uniform constant c > 0 such that (I6.12al) holds 
(resp., (16.12bl) holds). 

Proof. We only do the proof for (I6.12al) . the proof for (I6.12bl) is similar. Let 
Vh £ Hh,xn be a nonzero discrete field. Let <j>(vh) £ Hq (D) be the solution to the 
following Poisson problem: 


S7-(tV<j)(vh)) = V-(cu^), (/>{vh)\dD = 0. 

Note that this problem is well-posed owing to the assumed regularity and structure 
of £. Let us define v(vh) '■= Vh — Vcf>(vh)- This definition implies that 


V-(cv(v h )) = 0, Vx(w(w h )) = Vxv h , v(v h )xn\ dD = 0, 
so that v(vh) £ H xrl . We now bound ||tt/i||i, 2 (D) as follows: 


c\\ v h\\ 2 L 2 ( D ) < J (cv h )-v h dx = j (cv h )-(v h - Jh 0 v(v h ) + J£ 0 v(v h ))dx 

= / (cv h )-Jf 0 (v h - v(v h ))dx + / (i Ev h )-Jh 0 v(v h )dx 
J D J D 

= / (vv h )-Jh 0 V(<l>(vh)) dx+ / (cv h )-J£ 0 v(v h ) dx. 


J D J D 

Note here that we used that fff 0 Vh = v^. Then using the commuting property 
^{Jh t oH v h)) and since J^ 0 maps onto P ®(%), we infer that 


c\\v h \\ 2 L 2 (D) < / (ev h yy{Jl 0 <t>(v h ))dx+ / (<cv h )-Jf 0 v(v h )dx 

J D J D 

= / (cv h )-jf 0 v (v h )dx < c'\\v h \\ L 2 { D)\\JhO v ( v h)\\L^D)- 
J D 

The uniform boundedness of J7^ 0 on L 2 (D) and (|6.11a[) with s = 0 imply 

IKIIr 2 (r>) < c\\Jf 0 v(v h )\\ L 2 {D) < c'\\v(v h )\\ L 2 {D) < c ,, \\Vxv h \\ L 2 {D ). 

This concludes the proof. □ 
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